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Abstract 

A lattice field theory approach to the statistical mechanics of charged polymers in elec- 
trolyte solutions [S. Tsonchev, R. D. Coalson, and A. Duncan, Phys. Rev. E 60, 4257, (1999)] 
is applied to the study of a polymer chain contained in a spherical cavity but able to diffuse 
into a surrounding gel. The distribution of the polymer chain between the cavity and the 
gel is described by its partition coefficient, which is computed as a function of the number of 
monomers in the chain, the monomer charge, and the ion concentrations in the solution. 



1 



1 Introduction 



The problem of partitioning of a polymer chain confined to move within large cavities 
embedded into a hydrogel has received attention with some interesting experiments 
performed by Liu et al. [1]. These authors investigated the so called "entropic 
trapping" phenomenon, which describes the preferential localization of a polymer 
chain within large cavities embedded in a hydrogel due to the larger conformational 
entropy of the chain in them. Therefore, this phenomenon has been suggested as a 
basis for potential new methods of polymer separation. In this context, the prob- 
lem is also relevant and can lead to better understanding and possible improvement 
of many existing separation methods, such as membrane separation, filtration, gel 
electrophoresis, size exclusion chromatography, etc. [2], all of which utilize the de- 
pendence of macromolecular mobility through a network of random obstacles on 
molecular properties, such as molecular weight, monomer charge, electrolyte com- 
position, etc. The practical importance of these techniques has motivated a number 
of investigations of polymer separation between cavities of different size [3-7] . 

In our earlier work [6, 7] we used lattice field theory calculations to study polymer 
separation between two spheres of different size — a simplified model of the more 
complicated system of the polymer moving between large cavities embedded in a 
hydrogel, with the larger sphere playing the role of the cavities and the small sphere 
corresponding to the connecting channels in the gel. We investigated the dependence 
of the partition coefficient K, defined as the ratio of the average number of monomers 
in the two respective spheres, as a function of the total number of monomers in the 
chain, the excluded volume interaction between them, the monomer charge, and the 
concentration of electrolytes in the solution. Our results were qualitatively in accord 
with the experiments of Liu et al. [1] and with related computer simulations [8]. 

In this work we apply the lattice field theory approach to the more complex 
system of a polymer chain moving within a large spherical cavity embedded in a 
network of random obstacles. 
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In Section 2 of the paper, for continuity of the presentation, we review the lattice 
field theory of charged polymer chains in electrolyte solution [7]. In Section 3 we 
describe the Lanczos approach for finding the energy spectrum of the Schrodinger 
Hamiltonian problem (arising from the polymer part of the partition function [7]), 
and the resolvent approach for extracting the corresponding eigenvectors. Section 4 
describes the numerical procedure for solving the mean field equations of the system, 
and in Section 5 we present and discuss our results. In Section 6 we conclude our 
presentation. 
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2 Review of Lattice Field Theory of Charged Polymer Chains 
in Electrolyte Solution 



In Ref. [9] we derived the following functional integral expression for the full par- 
tition function of a charged polymer in an electrolyte solution with short-range 
monomer repulsion interactions 



Here, (3=1 /kT is the inverse temperature, e is the dielectric constant of the solution, 
e is the proton charge, A is a measure of the strength of the excluded volume inter- 
action, x an d uj are auxiliary fields, c±=e f5fl± / \± with fi± and A± being the chemical 
potentials and the thermal deBroglie wavelengths for the ions, respectively. The 
polymer part Zschrix^) i n (1) refers to a Euclidean-time (T=M=total number of 
monomers) amplitude for an equivalent Schrodinger problem based on the Hamilto- 
nian 



where a p is the Kuhn length and p is the charge per monomer. The mean-field 
equations corresponding to the purely-real saddle-point configuration fields Xc = % 
oj c = iuj are obtained by setting the variational derivative of the exponent in the full 
functional integral (1) to zero. For the case of a polymer with free ends (the only 
situation considered in this paper), the polymer amplitude -Zschr can be written in 
terms of sums over eigenstates of H as follows [7]: 



Z = f Dx{r)Duj(r)e^ J 



x Axdr-± Juj(f) 2 df+c + J e iel3x dr+c- J e- iel3x dr^ 



'schriX,^). (1) 



(2) 




n 



-ME, 



J2A 2 n e- ME " = e F v° l , 



(3) 



n 



where E n is the n-th energy eigenvalue, 




(4) 



4 



and 

F pol = \n(^Ale- ME ^ (5) 

is the negative of the polymer contribution to the free energy. Thus, the mean-field 
result for the negative of the total free energy is 

F = J dr^L |v Xc | 2 + K? c + c+eto*' + c_e~^} + F pol { Xc ^ c ) . (6) 

Varying the functional (6) with respect to the fields Xc, u c one obtains the mean-field 
equations 



e -l 



4ire 



V 2 Xc (f) = c+e^ Xe(f) - c-e-P***^ - pp(f) , (7) 



a 2 - 

-fV 2 V n (r) = Xp(f)^ n (f)+(3pe Xc (mn(r)-(E n -V m (f))^ n (f), (8) 



where p, defined as 

An^nAm^m ( p -ME n _ p -ME m \ 
(ff\= - ^ n ' m E "~ E ™ V > (Q) 

is the total monomer density. The equations presented here apply for polymer chains 
of arbitrary length, provided all (or a sufficient number) of states are included in 
the sums above. The single-particle potential V m (r) has been included to enforce an 
exclusion region for the monomers [9]. Note that the parameters c± are exponen- 
tials of the chemical potentials p± for positively and negatively charged ions. The 
numbers of these ions must be fixed by suitably adjusting c± to satisfy the relations 

n ± = c ± dl J {Z) =c ± [e^df. (10) 
oc± J 

The advantage of working with F is that, as shown in Ref. [7], it has a unique 
minimum, and thus, can be used to guide a numerical search for the mean elec- 
trostatic and monomer density fields. Once the mean fields have been computed, 
the defining relation \nZ=F( Xc , w c ) can be used to obtain free energies of various 
types. For example, the Helmholtz free energy A (corresponding to fixed numbers 
of monomers and impurity ions) is given by 

f3 A = n + lnc + + n_ In c_ — F( Xc , w c ) • (H) 



Following the procedure of Ref. [9], we now move from the continuum to a 
discrete 3-dimensional lattice by rescaling according to 

f(f) -> (3exc(r) 
^v(r) - af 2 ^ N {r) 

and multiplying Eq. (7) by af (aj being the lattice spacing). This leads to the 
following discretized version of equations (7) and (8) on a 3D lattice: 

a Anmfm = - 7_e~ / ™ - pp H , (12) 



m 



^X)^™*^™ = —^-Pn^N,n+Pfn^N,n-E N ^ NtH , (13) 



oaf - a; 



where 



* = (15) 

and the wavefunctions are dimensionless and normalized according to 

E^v,n = i; (16) 

n 

thus, the density p^j sums to the total number of monomers, M. 
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3 Extraction of Eigenspectrum and Eigenfunctions for Poly- 
mer Effective Hamiltonian 

The simultaneous relaxation solution of Equations (12) and (13) requires a rapid and 
efficient extraction of the eigenvalues and low-lying eigenvectors of the operator H, 
which amounts — once the problem has been set up on a discrete finite 3-dimensional 
lattice — to a large sparse real symmetric matrix. 

We have found it convenient to use distinct algorithms to extract the low-lying 
spectrum and eigenvectors of H (typically we need on the order of 10-30 of the 
lowest states for the shortest polymer chains studied here, while for the longest 
polymer chains only one to three states suffice). The eigenvalues are extracted using 
the Lanczos technique [10]. Starting from a random initial vector wo = vi, one 
generates a series of orthonormal vectors v±, i>2, ... by the following recursion: 

v n+1 = w n /p n , 

n — > n + 1 , 

a n = (v n , Hv n ) , 

w n = (H - a n I)v n - /3 n _iv n _i , 

Pn = \J(w n ,W n ) , 

where a n , f3 n are real numbers, with fi = 1 and vq = 0. The matrix of H in the basis 
spanned by v n is tridiagonal with the number a n {(5 n ) on the diagonal (respectively, 
super/sub diagonal). Carrying the Lanczos recursion to order N, diagonalization of 
the resulting NxN tridiagonal matrix leads, for large N, to increasingly accurate 
approximants to the exact eigenvalues of H. The presence of spurious eigenvalues 
(which must be removed by the sieve method of Cullum and Willoughby [11]) means 
that typically a few hundred Lanczos steps must be performed to extract the lowest 
30 or 40 eigenvalues of H (for dimensions of H of order 10 5 as studied here) to 
double precision. 

Once the low-lying spectrum of H has been extracted by the Lanczos procedure, 
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as outlined above, the corresponding eigenvectors are best obtained by a resolvent 
procedure. Supposing \ n to be the exact nth eigenvalue of H (as obtained by the 
Lanczos method), and ip n the corresponding eigenvector, then for any random vector 
■0 ran with nonzero overlap with ip n , the vector obtained by applying the resolvent 

FT^ran (17) 

A n + e — ti 

is an increasingly accurate (unnormalized) approximant to the exact eigenvector ip n 
as the shift e is taken to zero. A convenient algorithm for performing the desired 
inverse is the biconjugate gradient method (see routine linbcg in [12]). We have 
found the combination of Lanczos and conjugate gradient techniques to be a rapid 
and efficient approach to the extraction of the needed low-lying spectrum. 
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4 Solving the Mean-Field Equations for a Polymer Chain 
Confined to Move within a Spherical Cavity Embedded 
in a Gel 

Equations (12) and (13) are solved simultaneously using the following relaxation 
procedure [9]. First, the Schrodinger Eq. (13) is solved for fn—0 and ignoring the 
nonlinear (monomer repulsion) potential term. The resulting Stan's and correspond- 
ing energy levels (wavefunctions and energy eigenvalues of a particle confined to 
the cavity in a gel system) are used to calculate p^, then the Poisson-Boltzmann Eq. 

(12) is solved at each lattice point using a simple line minimization procedure [13]. 
The process is repeated and the coefficients r y± are updated after a few iterations 
until a predetermined accuracy is achieved. Then the resulting is used in Eq. 

(13) , which is solved using the Lanczos method [10] for a new set of ^n,hs to be 
used in calculating an updated version of the monomer density p^. This density is 
then inserted into Eq. (12) and a new version of is computed. For numerical 
stability, the updated inserted into Eq. (13) is obtained by adding a small frac- 
tion of the new fn (just obtained from Eq. (12)) to the old one (saved from the 
previous iteration). The same "slow charging" procedure is used for updating in 
the nonlinear potential term of the Schrodinger equation (13). 

This numerical procedure has been applied to the system of a polymer chain 
moving within a cavity embedded in a network of random obstacles. We carve a 
spherical cavity of radius 10a; in the middle a cube with a side-length of 40cfy on 
40 3 lattice. The Kuhn length, a p = 2a«, and in absolute units a p = 5A. Then the 
random obstacles are created by randomly selecting 20% of the remaining lattice 
points in the cube outside the carved sphere to be off limits for the polymer chain. 
Thus, the random obstacles take 20% of the gel, that is, 80% of the gel volume 
plus the cavity volume is available for the chain to move in. On the other hand, 
the impurity ions are free to move within the whole volume of the system. The 
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monomer repulsion parameter A is fixed throughout the computations through the 
dimensionless parameter ( by the following relation: 

C = 4v4, (18) 

and the dimensionless parameter £ = 5. 
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Figure 1 : In K vs M for varying monomer charge p and fixed number of negative impurity ions 
n = 600, which corresponds to molar concentration C w 0.996M. 

5 Numerical Results and Discussion 

We have computed the log of the partition coefficient K = (Mi) / (M2), where (Mi) 
and (M 2 ) are the number of monomers in the spherical cavity and the remaining 
gel, respectively, as a function of the total number of monomers in the system, 
M = (Mi) + (M2), for varying monomer charge p and varying number of ions in the 
system. In Fig. 1 we show the plot of In K vs M for two different monomer charges, 
p = —0.1 and p = —0.2 (all in units of e), and a fixed number of 600 negative 
impurity coions in the system, while the number of the positive counterions is fixed 
according to the condition for electroneutrality. We see that the partition coefficient 
K increases with M for only the shortest polymer chains, goes through a turnover, 
and from then on decreases continuously as the number of monomers is increased. 
As in our previous work [6, 7], we observe that smaller monomer charge leads to 
higher partition coefficient, due to the weaker repulsion between the monomers. In 
Fig. 2 we show how the partition coefficient varies as we vary the number of negative 
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Figure 2: In if vs M for varying number of negative impurity ions n and fixed monomer charge 
p = -0.1. 

impurity ions in the system. As expected, the higher number of ions leads to better 
screening of the monomer charges, hence less repulsion and larger K. Qualitatively, 
this behavior is similar to what we observed in our previous work [6, 7]; however, 
the partition coefficient shown here decreases for almost the whole range of M, and 
is, in fact, much larger than the coefficients reported earlier [6, 7] for partitioning 
of a polymer chain between two spheres. This can be explained as a result of the 
much smaller voids that arise between the random obstacles in the gel outside of the 
spherical cavity, compared to the smaller of the two spheres treated in [6, 7], which, 
in the Schrodinger language, means that, even though the volume available for the 
polymer chain outside of the spherical cavity is much greater than the volume of the 
cavity itself, the energy levels of the excited states which lead to a higher monomer 
density outside of the cavity are too high (due to the strong confinement in the 
narrow voids), so that the chain is largely confined to the cavity. Only for the cases 
of very large M do we observe a non-negligible monomer density outside of the cavity. 
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Figure 3: The average radial density p(r) as a function of the distance from the center of the 
spherical cavity r for the three sets of parameters considered here in the case of M = 40. 

This is illustrated in Figs. 3 and 4, where we plot the averaged radial density of 
monomers starting from the center of the spherical cavity for the three different 
sets of monomer charge and impurity ion concentration parameters presented here. 
In Fig. 3 we plot the radial density for the case of relatively small number of 
monomers, M = 40, and we see that virtually all of the monomers are confined 
to the spherical cavity, while in Fig. 4, which represents the case of M = 300, we 
observe a small but non-negligible contribution to the monomer density from the 
region outside of the cavity. In Figs. 5 and 6 we show the plots of the electric 
potential f(r) corresponding to the parameters of Figs. 3 and 4, respectively. We 
can qualitatively compare the results from Figs. 3-6 to our previous results in [9], 
where we computed the monomer density and the electric potential for a charged 
polymer chain confined to move within a sphere. It is clear that in both cases 
the shape of the monomer density distribution and the electric potential are quite 
similar, which is an illustration of the fact that the spherical cavity embedded in 
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Figure 4: Same as in Fig. 3, but for M = 300. 



the gel does indeed act as an "entropic trap" for the polymer chain, and for most 
of the range of reasonable physical parameters the system behaves approximately 
as a polymer chain in a spherical cavity. In Figs. 5 and 6 we see that, for the case 
of lower counterion numbers, the potential f(r) drops to negative values at large 
radial distance. Nevertheless, it does approach (up to finite lattice size corrections) 
zero slope, or equivalently, zero electric field, consistent with the overall electrical 
neutrality of the system. 
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Figure 5: Electric potential /(r) for the parameters corresponding to Fig. 3. 
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Figure 6: Electric potential /(r) for the parameters corresponding to Fig. 4. 
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6 Conclusions 



We have applied a previously developed lattice field theory approach to the statistical 
mechanics of a charged polymer chain in electrolyte solution [9, 6, 7] to the problem 
of a charged polymer chain moving in a spherical cavity embedded in a gel. This 
problem is more relevant to real experimental situations involving charged polymer 
chains in a complex environment than the two-sphere problem studied by us earlier 
[6, 7]. The results of this work demonstrate the capability of the approach to treat 
more complex systems of arbitrary shape in three dimensions, and also confirm the 
expectations that a large spherical void carved out from a network of random obsta- 
cles can act as a "trap" for polymer chains, and therefore, may serve as a prototype 
for new methods of polymer separation based on macro molecular weight, monomer 
charge, and/or electrolyte composition. The results presented here confirm our pre- 
vious contention [9, 6, 7] that chains with smaller monomer charge would be easier 
to separate by a technique exploiting the idea of "entropic trapping." Similarly, for 
chains with fixed monomer charge, a better separation would be achieved in solu- 
tions with higher impurity ion concentration — a parameter which is typically varied 
in the laboratory. 

It is important to note that the method used here is based on the mean field 
approximation, and therefore, the results should be considered only as qualitative. 
Nevertheless, one can expect that the long range of the electrostatic interaction and 
the strong confinement of the polymer chain inside the spherical cavity would result 
in weakly fluctuating density and electrostatic fields and would make the mean field 
approximation reliable [14]. 
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